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I. INTRODUCTION 



The large N c limit, where N c is the number of colors, is a useful device to understand 
many systematic features of baryon properties |],|2|], such as the l/N c scaling of various 
physical quantities. In a series of papers, Dashen and Manohar |||4]], and Jenkins || have 
discussed the 1/N C structure of baryon properties, and the framework for combining chiral 
symmetry with the large N c aspects of QCD has been developed by many authors @-12]. In 
Refs. |f7|,|T0l|, the baryon octet and decuplet mass spectra were discussed in this framework 
and the baryon mass relations were derived. However, although those works successfully 
reproduce mass relations at tree level, they do not compute all possible terms allowed by 
the chiral and large N c expansions. 

The baryon mass spectrum was re-examined in conventional baryon chiral perturbation 
theory by Borasoy and Meissner ||T3l . |T4| . To compute the baryon masses to order m 2 q , where 
m q is the quark mass, the decuplet degrees of freedom are integrated out to give counter 
terms, and some low-energy constants are determined from resonance saturation. However, 
when we work with the 1/N C expansion, the octet and decuplet states are degenerate at the 
leading order, and the decuplet fields must be treated explicitly. 

In this paper, we re-examine the baryon masses in chiral perturbation theory taking into 
account the l/N c counting based on the techniques developed in the literature, e.g., in Refs. 
P Pf. This enables us to investigate the 1/N C structure of the baryon properties and the 
meson-baryon interactions in a systematic way. The baryon axial current matrix elements 
and the strangeness contribution to the nucleon mass are computed as well. Some of these 
topics were discussed in the literature |^, [10 , 11| focusing on the leading order terms in 1/N C 
expansion (up to one loop corrections). In this paper, we perform the calculations to the 
next orders and discuss a difficulty which arises in computing the one loop corrections in a 
way which is consistent with the 1/N e expansion. This paper is organized as follows. In the 
next Section, we briefly discuss the formalism of this approach. In Section III, we compute 
the baryon axial current up to one-loop corrections. The baryon mass formulas are given 
in Section IV. The one-loop corrections to the baryon masses are calculated in Section V. 
We discuss the strangeness contribution to the nucleon mass and the sigma term in Section 
VI and then finish with a summary and conclusion in Section VII. Explicit expressions of 
baryon wave functions and some detailed formulas are given in Appendices. 



II. FORMALISM 

We start with a brief review of the construction of baryon states in the large N c limit, 
referring to Refs. 0,[n],[nj for further details. The baryon states with N c quarks can be 
written as follows: 

I R\ — U a l a \--- a N c ^ Nc A-i...A N J f \n.\ Icy -i \ 

in particle number space, where aj are the flavor indices, q.{ the spin indices and Ai the 
color indices. The quark creation and annihilation operators and q satisfy the usual anti- 
commutator relations for fermions. The symmetric tensor B is characteristic of each given 
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baryon wave function. Since the baryons are in color-singlet states, however, it is more 
convenient to work with 



\B)=B^-- a ^al iai ...al NcaNc \0), (2.2) 

by dropping the explicit color indices, where the operators and a are bosonic operators 
satisfying the usual commutator relations. For short-hand notation, we label the quark 
operators as 

a 4 = ai d i, a 5 = a sh a 6 = a sl , (2.3) 

so that oc[ creates w-quark with spin-up, and so forth. 

There is an ambiguity when we extrapolate the physical baryon states to large N c . As 
in the literature, we keep the spin, isospin and strangeness of baryons as 0(1) in the large 
iV c limit. For example, the nucleon state in large N c limit has spin 1/2, isospin 1/2 and no 
strangeness. This can be done by acting with spin-flavor singlet operators on the physical 
baryon states. For example, the proton spin-up state can be written as 

\p,+l)=C N a\(Al) n \0), (2.4) 

where n = (N c — l)/2 and Cn is the normalization constant. The spin-isospin singlet 
operator A\ is defined as 

A\ = a\a\ - a\a\. (2.5) 

One can easily verify that this state reduces to the usual quark model state in the real world 
with N c = 3. The complete list of the baryon octet and decuplet states can be found in 
Appendix A. 

Next we define a one-body operator {X} in spin-flavor space as 

{X} = aLXJV (2-6) 

so that its expectation value on baryon states is at most of 0(N C ). In a similar way, one 
can define 2-body operators and 3-body operators {X}{Y}{Z}, and so on. Then, 

it is found that the coefficient of an r-body operator is at most 0(1/N^ where £ is 
the number of inner quark loops P,|T6|. This enables us to treat the coupling constants as 



0(1) quantities in the large N c expansion by writing the iV c -dependence of the operators 
explicitly. 

By direct evaluation, one can see the well-known commutator relation, 

\{X} 1 {Y}] = {[X,Y]}. (2.7) 

Note that the left-hand side is naively a two body operator whose expectation values can 
be of O(N^), whereas the right-hand side is a one-body operator whose expectation values 
are of 0(N C ) at most. This means that the order of an operator in 1/N C counting reduces 
when we have a commutator structure as in Eq. ( |2.7|) . This plays an important role in the 
large N c analyses of the baryon properties. 

We will discuss the explicit forms of some operators which appear in the calculation of 
baryon axial currents and masses in the next Sections. 
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III. BARYON AXIAL CURRENTS 



A. Tree level 

Our starting point is the chiral meson-baryon effective Lagrangian. Baryon matrix ele- 
ments of this Lagrangian involve the meson-baryon interaction in the following form: 



(4*> = 9 (B\{A^}\B) + A (B\{A»}{a„}\B) + ..., (3.1) 



where is the baryon spin matrix,^] and the dots denote higher order terms. The axial field 
is defined as 

A,= l -(td^-t;%t), (3.2) 

where £ = exp(ill//) with the meson decay constant /. The SU(3) matrix field II represents 
the octet of pseudoscalar Goldstone bosons. It is defined as 

IT = ^\ a n a , (3.3) 



with the usual Gell-Mann matrices A a (a = 1, . . . , 8). In Eq. (|3.1| ), the N c factors of operators 
are given explicitly, and the coupling constants g and h are of 0(1) in the 1/N C counting. 
Then the baryon axial current J| „ reads 



Jl, = f {T a ^} + ^ {f a }K>, (3.4) 



from the Lagrangian Q3.1p with 



f a = \(Z\ a ? + ?\ a Z). (3-5) 
This gives its matrix elements as 

(B'\J^\B) = afaBUB'MuB. (3.6) 
where ub is the Dirac spinor of the baryon and 

°& B = 9(B'\{±\ a cr 3 }\B) + ^-(B'\{l\ a }{a 3 }\B) } (3.7) 

at the tree level. 

By using the wave functions given in Appendix A, we can compute the baryon axial 
current straightforwardly. A naive investigation of each term gives that, despite the 1/N C 



1 In the baryon rest frame, = (0, a) with the usual Pauli matrices a 1 . 



4 



factor, the h term contribution is expected to be of the same order as that coming from the 
g term. This is because the h term contains a 2-body operator whose expectation value can 
be 0(7V 2 ), thus the leading order of the h term contribution can be of N c . However, close 
inspection shows that the g term contribution dominates, because the h term contains the 
operator {<r M } and our baryon wave functions satisfy {cr^} ~ 0(1). 
The explicit forms of the relevant operators are 

{|A 1+i2 o- 3 } = a\a 3 - c4a 4 , { |A 4+45 (X 3 } = a\a 5 - a\a & , 
{|A 1+l2 }{a 3 } = a{a 3 + «la 4) UA 4+i5 }{a 3 } = a\a 5 + a\a 6 , (3.8) 

so that we obtain 



< + ' 2 = i^+2 ) + A 

1-N2 _ ^,l+«2 _ 9 



«as 2 = 4 + +a = ;rW(*e - W + 3), 



a 



7(1- 



9 TV, 



c 



V2h 



and 



4^- = = ^(N c + 1) + (3.0) 



otf B = -|V JV '« + 3 -2^V-V.- + 3. 



A= 2^3 2iV c v 



These results show that the h term contributions are suppressed as compared to those of 
the g terms as we discussed above. We can also find that the leading order of a 1 ^? is 
0(N C ), whereas a^f, which changes the baryon strangeness, is 0(y/N^). This shows that 
the strangeness-changing (i.e., AS* ^ 0) baryon axial currents are suppressed as compared 
to the strangeness-conserving (i.e., AS = 0) baryon axial currents by 0(y/N2). This can 
be understood from Eq. (|3.8j ) by noting that the number of u,d quarks in the baryon wave 
functions is 0(N C ) whereas that of s quark, i.e., strangeness, is 0(N®). For example, in the 
case of «p+* 2 , acting with a 3 (or a 4 ) on the baryon state gives the factor N c , and the inner 
product of initial and final baryon wave functions with the proper normalization constants 
gives 0(1), so that ot^f 1 is 0(N C ). However, for a p j^ 5 , the action with a$ (or ag) gives 0(N®) 
because our baryon wave functions have the strangeness of 0(N®). Since the normalization 
constants of nucleon and A are 0(1/N C ) and 0(l/y/N^), respectively, we have an additional 
factor y/Nc in the calculation of at/C 5 , which implies that the order of a 4 .^ 5 i s O(v^Vc). 
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Since the contributions of the h term are suppressed as compared to those of the g term 
by 0(1/Nc) for the AS = axial currents and by 0(1/N C ) for the AS = 1 axial currents, 
we can neglect the h term up to next to leading order. At this order, when we fix N c = 3, 



we can recover the quark model relation [10 



D = g, F =ld, (3.H) 



by comparing with the results of the baryon chiral perturbation theory [|17],|]1| i n addition 
to the quark model predictions 

C = -2g, H = -3g, (3.12) 

for the octet-decuplet-meson and decuplet-decuplet-meson coupling constants, C and TC, 
defined as in Ref. [|l£j. When we go further in the 1/N C expansion, we must include the h 
term, and we get the modified relations, 

2q + h , 
D = g, F = JL ^-> (3-13) 

as found in Ref. [TO. 

We also compute a B , B by using 

{§AV} = A=(#i - N 2 + N 3 - iV 4 - 2N 5 + 2N 6 ), 

{lA 8 }{or 3 } = A=(iVi + N 2 + N 3 + N 4 - 2iV 5 - 2JV 6 ), (3.14) 

where we have introduced Nt = aqai. This leads to 

H q h 



VP 



2v/3 2\/3' 

g_ h ( 3 
~v / 3 + 2 v / 3V N c 
8 9 h (i 3 



2 2^ V iVj V ; 

From these results we find that the leading order of a%, B is 0(A r °) and that the h term 
provides a leading contribution together with the g term. This is because the expectation 
values of {f A 8 }{<7 3 } are 0(N C ) whereas those of {|A 8 a 3 } are 0(1). So we conclude that in 
order to get a consistent result on a%, B , one should consider n-body {n > 3) operators in 
general unless their coupling constants are suppressed. From the fitted values of D and F, 
one can estimate g = 0.61 ~ 0.8 together with h = —0.02 ~ —0.1, which shows that h is 



indeed small, less than 15% of g, but with opposite sign fllQ| . Therefore, one should keep in 
mind the contributions from n-body {n > 3) operators in the calculation of a B , B . We have 
a similar situation when we compute the 77-meson loop corrections to the baryon masses in 
Section V. 
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B. One-Loop Corrections 



The one-loop corrections to the baryon axial current in large N c chiral perturbation the- 
ory as shown in Fig. [I] have been discussed in Refs. 0,1,0- Naively, these loop corrections as 
they stand are not consistent with the 1/N C expansion. From the meson-baryon interactions 
([3.1D, each vertex is related to an operator X ta defined as 



X™ = <7{±AV} + A { i A «}{^} ; (3.16) 



c 



with spin index % and flavor index a. This shows that the meson-baryon coupling is of order 
N c . Because of the presence of 1/f which scales as 1/yiVc, each vertex has a factor \TN~ C . 
Then from Fig. |l|(a), it is evident that the one- loop correction is O(N^) when olb'b is of 
order N c . Thus the loop correction dominates the tree level. However, we have to consider 
the wave function renormalization terms (Fig. ||) in the loop calculation. When combined 
with Fig. [I], this gives the commutator structure to the baryon axial current operators, 
which implies that the actual order of the one-loop correction is 0(N®) when chb>b is 0(N C ). 
This suppression was proved up to 2-loop order in Ref. [^(J which concludes that the 2-loop 
corrections are suppressed by 1/N% as compared to the tree level values. 

Explicitly, the one-loop correction to the baryon axial current from Fig. |TJ(a) is given by 
the following expression: 

where vn^i{= m n ,mK,rn v ) is the meson mass in the loop. When combined with the wave 
function renormalization factor Zb from Fig. |2|, 

this gives the renormalized baryon axial current operator in the form of 

X ia +^I%[X^>,[X ia ,X vb % (3.19) 

where 

j oh , f d k 1 ki^ki/ , * 

= "V {2TxY{k-vfml b -k^ (3,20) 

Finally, we get the one-loop correction to the baryon axial current operator as 

2 / 2 o \ 2 2 

ryw m W f i m bb' | I T Vjb fyia vjb'll m bb' i_ m bb' r ^i,bb' (rtn-\„\ 

6X = 3Wl ^ + sJ [X AX ' X ]] ~3W ^ ' ( } 

by evaluating the loop integral using dimensional regularization with the regularization scale 



A (see, e.g., Ref. Plj.), where 
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Of = 9{[^\ [|A b ', §A°K} + -{[±A fe , [|A b ', |A°]]}K}, (3.21b) 

which comes from Fig. |l](b). So the one-loop correction to the baryon axial current matrix 
elements are obtained as 

2 2 2 

5a B , B = l3 BIB ^^\n— + p B , B —^, (3.22) 

where II stands for tt, K and rj mesons. 

The explicit results of flgjg and ^ B ^ B with g terms are given in Appendix B. From these 
results, we can see that the one- loop corrections are 0(1/N C ) at most since f 2 scales like N c . 
Furthermore, the corrections from Fig. |](b) are of the same order as those of Fig. |](a). 

IV. BARYON MASSES AT TREE LEVEL 

In this Section we investigate the baryon masses at tree level. To estimate the baryon 
masses simultaneously in the 1/N C expansion and the chiral expansion, we must specify the 
relation between l/N c and the pseudo-Goldstone boson mass m n . In this paper, we use 
muSM ~ 0(1) where 5M is the octet-decuplet mass difference. This gives m n ~ 0(e) 
and 1/N C ~ 0(e), where e stands for a small expansion parameter. 2 A priori, there is 



no constraint on the relationship between mn and N c . In fact, the authors of Ref. fll 
used m^SM ~ 0(1). However as we shall see below, the leading order correction to the 
degenerate baryon mass in the large N c limit is proportional to m^Nc and we count it as 
0(e). This is consistent, given that mn ~ 0(e 2 ) in accordance with the chiral expansion, 
and N c ~ 0(e^ 1 ). We will compare our results with those of Ref. [11 before calculating the 
one-loop corrections. 

The matrix elements of the effective Lagrangian which contribute to the baryon mass 
can be written as 

(C B )=j:(B\C®\B), (4.1) 

i 

where C^g represents that part of the Lagrangian which can give a contribution of 0(e l ). 
Explicitly, these terms are 

4 ff 1} = -oo{l}, (4.2a) 

4S = -^HH-&iW, (4.2b) 

4S = -^tr(m){l}, (4.2c) 



2 This is consistent with the expansion of Ref. 23], where the A-nucleon mass difference is treated 
as small parameter with the pion mass. 
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£2 = -h.{moi}{oS} - Cl {m 2 } - ^ {m}{m}, (4.2d) 

= -^tr(m)K}{^} - ^tr(m 2 ){l}, (4.2e) 
£g = ~{mV}{ai} - ^H^} - ^W{^}{^} 

- ^{m 3 } - -|{m 2 }{m} - A {m}{m}{m} , (4 . 2 f) 

up to 0(e 5 ), where 

m = B @m q g + Zm g £). (4.3) 

We make use of the standard relations between B and squared pion and kaon masses, 
m 2 = 2B m and m 2 K = B {rh + m s ), where m is the average mass of u and d quarks and 
m s the s-quark mass. The quark mass matrix m q is given by 

m q = mlA + m s S, (4.4) 

where 

W = diag(l,l,0), S = diag(0,0, 1). (4.5) 

Throughout this work, we assume SU(2) isospin symmetry with m u = m^. Then, up to 
0(e 5 ), there are 15 low energy constants that should be determined from experiments. 
However, one can find that 6 terms give identical contributions to all baryon masses so 
that 9 parameters remain which determine all baryon mass differences. In the following, we 
discuss the baryon masses at each order of e. 

From the Lagrangian ( |4.2a| ), all octet and decuplet baryon masses are degenerate at 
leading order, which gives the baryon mass operator, 

Mf = a N c , (4.6) 

where the parameter ao sets the scale as a "mass per color degree of freedom" . 
To the next order, the correction to the mass formula reads 

5mS ] = ^-{<T j }{cr j } + 2B mN c h. (4.7) 

The a\ term gives the splitting between octet and decuplet while all states within the 
octet and decuplet are still degenerate. Although the original form of Eq. (|4.2b[) includes 
the operator {S}, the resulting baryon masses do not depend on strangeness since the 
expectation values of {S} for our baryon states are of 0(1) so that its contribution appears 
at the next higher order. Thus, at O^e 1 ), we get 

5M^ = J- ai + 2B mN c b 1 , 

8M$ = —ai + 2B mN c b u (4.8) 
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where M 8 and Mi denote the baryon octet and decuplet masses, respectively. 



At 0{e 2 ) there are two contributions. One is from L 2 S of ( f4.2c|) and the other is from 
the remaining part of the b\ term of Cl s : 

5M { z ] = 2B (m s - m)6i{5} + 2B (m s + 2m) ai . (4.9) 

It is clear that the a\ term gives the same mass shift to all baryons. The dependence of the 
h\ term on strangeness results from the SU(3) flavor symmetry breaking and vanishes in the 
flavor SU(3) limit. Therefore, up to this order, the baryon mass depends on total spin and 
strangeness, but the A and the £ are still degenerate. 
The mass corrections at 0(e 3 ) can be obtained as 

SMP = ^mb 2 {a j }{a j } + ^(m s - m)b 2 {S^}{^} + N c (2B m) 2 ( Cl + c 2 ). (4.10) 

The b 2 term involves two operators. One of them depends on the total baryon spin and the 
other depends on the spin of the strange quark(s). As a result, the S decouples from the A 
at this order. Up to this order, we have 4 operators in the baryon mass formula, namely, 
{1}, {cr^Hcx- 7 }, {<S} and {Sa^}{a^}. The c\ and c 2 terms give the same contributions to all 
baryons. The matrix elements of the operators can be evaluated using 

{S} = (N 5 + N 6 ), 
{S}{S} = (N 5 + N 6 ) 2 , 
{Saj}^} = 2(N 5 + N 6 ) + {N 5 - iV 6 )[(iV 1 - N 2 ) + {N 3 - N 4 ) + {N 5 - N 6 )] + 4iV 5 iV 6 
+ 2{aia 2 a\ctQ + a\a 2 a^a\ + a 3 a\a\a§ + a\a 4 a$a\), 
{^{aj} = 2[(Ni + N 2 ) + (N 3 + N 4 ) + (N 5 + N 6 )\ + 4(NiN 2 + N 3 N 4 + N 5 N 6 ) 
+ [(Ni - N 2 ) + (N 3 - N 4 ) + (N 5 - N 6 )] 2 

+ A(a\a 2 a 3 a\ + a\a 2 a^,a\ + aia^a^a^ + aia 2 a\a 6 + a 3 a\a\a§ 
+ a 3 a\ala 6 ), 
{Sa }{Ua } = [(Ni - N 2 ) + (N 3 - N 4 )](N 5 - N 6 ) 

+ 2(ala 2 a^al + aia 2 a\a§ + a^a^a^al, + a 3 a\a\a§). (4-H) 

All the matrix elements needed to calculate the baryon masses are given in Table |. 
The explicit expression of mass corrections at 0(e A ) reads 

5M^ = (2B ) 2 [(m 2 s - m 2 ) Cl + 2m(m s - m)c 2 ]{S} 

+ ^2B (2m + m s ){a j }{a j } + fr{2B Q ) 2 (2m 2 + m 2 ). (4.12) 

N c 

The combination of c\ and c 2 terms depends on the strangeness, and the a 2 term gives the 
next order contribution to the decuplet-octet splitting. Therefore, all of the above terms 
can be absorbed into the formulas valid up to 0(e 3 ). 
Then, up to this order, we have three mass relations, 
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(Ml) M A - M N Ms-, - M s • ~(M S - M A ), 

(M2) 3 Ma + M E - 2M V - 2M H = 0, 

(MS) Mo - M=* = M H * - M 2 * = M 2 * - M A , 



(4.13) 



where (Mi) is the hyperfine splitting rule, (M2) the Gell-Mann-Okubo (GMO) relation and 
(MS) the decuplet equal spacing (DES) rule. 

The 0(e 5 ) correction to the baryon mass has a more complicated form: 



5M 



(5) 



N, 



(2B ) 2 (m s -m) 2 c 2 {S}{S} 



c 3 



+ -f(2B ) 2 [m 2 {a j }{a j } + (m 2 s - m 2 ){Sa j }{a'} 
+ ^(25 ) 2 [mVH^} + {m 2 s - m 2 ){Sa^} 
+ ^(2S ) 2 [m 2 + m(m s - m){Sa j }{a j } 



m) 2 {^'}{5(7 J '} 



+ iy c (d 1 + d 2 + d 3 )(2B m) 3 . (4.14) 

There are more terms including C5 and dx,2,3> but they give contributions only at higher 
orders. The mass formula ( 4.14 ) includes the operators and {Wo"- ? }{ l Scr J } in addition 

to the operators that appeared already at the lower order. Because of these new operators, 
the mass relations (M2) and (MS) of Eq. ( |4.13j ) are modified, whereas (Mi) is still valid. 
Instead of (M2) and (MS), we find improved GMO and DES rules [|,[K]]: 



(M2) 3M A + M s - 2(Mv + M H ) = (M E * - M A ) - (M n - Ms*), 
(M^) (M n - M H *) - (Ms. - M s .) = (M H . - M s .) - (M s . - M A ) 



(4.15) 



which work better than (M2) and (MS). Empirically, the left and right hand sides of (Ml) 
give (293 = 308), and (M2) and (MS), respectively, lead to (27 = 0) and (139 = 149 = 152), 
whereas (M2 ! ) gives (27 = 11) and (MS') gives (—3 = —8), where the numbers are given in 
MeV. Combining these relations with (Ml ) gives 



Me 



- Ms = M s . 
(215 = 192) 



My 



(4.16) 



where the numbers show again the experimental values. Note that this is not an independent 
mass relation. The modified DES rule (MS') was first derived by Okubo \T2\ in the form of 

M n - M A = 3(M H * - M E .) , (4. 17) 

which is just a re-combination of (Mi), (M2') and (MS'). 

Since there are 6 different types of operators up to 0(e 5 ), we can write the mass formula 
in a compact form as 



M B = a + b{a j }{a j } + c{S} + d{Sa j }{a j } + e{S}{S} + f{Ua J }{Sa 3 }, 



(4.18) 
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where the c term comes in at 0(e 2 ), the d term at 0{e 3 ), and the e and / terms at 0{e 5 ). 
The best x 2 fits to the baryon masses up to 0(e 5 ) are shown in Table || and Fig. The 
best fit up to 0{e A ) is the same as that of 0{e 3 ). This is because the mass formula of 0(e 4 ) 
does not introduce any new operator. A reasonable baryon mass spectrum is already found 
at 0{e 3 ). Corrections from the 0{e 5 ) operators are evidently not so important. Note also 
that the coefficients of the operators involving S include a factor (m s — m) so that the c, d, 
e and / terms of Eq. (|4.18[ ) vanish in the limit of exact SU(3) flavor symmetry. 

Before proceeding to the loop corrections, let us compare our results with those of Ref. 
TT|, which uses different counting so that O(e') = 0{m q ) = 0(1/N C ). At the leading order, 



the authors of Ref. [nj obtained 5 mass relations: 



(M H - 


-M E ) 


-(M s - 


M N ) + 


^(Me-M a ) 


(Mn< 


-M s . 


) - (M 2 * 


— Ma) 


= (=-8) 


(M n - 


-M B *) 


- {Mb* 


- Ms*) 


= (=-3) 


(M s - 


-M N ) 


- (M A - 


M w ) = 


(=77), 




— Ma) 


-(M A - 


- Miv) = 


= (= -24) 



-13.5), 



(4.19) 

where the numbers in parenthesis on the right hand side are the empirical ones in MeV. The 
first three relations are re-combinations of {Ml), {M2) and {MS) and they are reasonably 
consistent with experiments. However, the deviations of the last two relations are larger 
compared with the first three relations. In our scheme, this discrepancy can be understood 
easily because the first three relations hold up to 0{e i ) and 0{e 4 ) whereas the last two hold 
only up to 0{e 2 ). 



V. LOOP CORRECTIONS TO BARYON MASSES 

The one-loop corrections to the baryon masses are obtained from the diagrams shown 
in Fig. |j. First, let us consider the diagram of Fig. |](a) without mass insertions to the 
intermediate baryon states, which corresponds to Fig. ||(a). At first glance, this one- loop 
correction appears to be inconsistent with the 1/N C expansion. Since each vertex carries a 
factor y/~N~ c , the one- loop correction is 0{N C ). A similar feature occurs in the case of the 
baryon axial current, where the wave function renormalization part must be included to give 
the proper commutator structure which is essential to be consistent with the 1/N C expansion. 
In the case of the baryon self energy, however, there is no other term that can lead to this 
commutator structure. Thus the one-loop correction is not suppressed as compared to the 
tree level baryon masses 0. In fact, this one-loop correction starts from 0{N C ), but we can 
see that the corrections of this order are the same to all baryons so that it can be absorbed 
in the term of the baryon mass. 

The one-loop baryon self energy is obtained as 

E B (.) = ±p\X-X~\B) J ^J^ ' mL (-*»*»), (5.D 
where u = p ■ v and X^ a is defined in Eq. ( |3.16|) . After evaluating the loop integral we find: 
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5M F 



nr 



n //on 



16tt/ s 



with 



(5.2a) 



(O u (B)) = (B\O n \B) = -(B^X^X^B), 

3 „ 



(5.2b) 



where a = 1,2,3 for the pion loop, a = 4, ... ,7 for the kaon loop, and a 
loop. The operator O n can be computed straightforwardly to give 



8 for the eta 



0* = f 



K = g< 



+ 



+ 



+ 



0"=g l 



yy2 / 8 \ 1 1 2 

-f + 2iV c - [N c + -j {S} + -{S}{S} - ~{v j }W j } + 3 {Soi}{oi} 
gh 



+ 3N [Nc ~ {S} + 2] i^>^> - {^H^} 



6JV? 



K}{a<} - 2{S<r i }{<r i } + 2{5} + {^}{^} 



^ [({5} + 1)M V} + (N c - 2{S} + 1){<S^}M 



3N? 



N c + (N c - 1){S} + {Sa*}{a 1 } - 2{S}{S} {a*} {a*} 



{S} + - -{Sa^Hai} + ±{o*}{oS} 

gh 



+ ^[Nc - S{S}} [{o>}{o>} - 3{So>}{o>} 



+ 



h 2 



18W 



(5.3a) 



(5.3b) 



(5.3c) 



There are several remarks concerning this result. As we discussed before, the pion loop 
correction O n includes the factor N 2 , which gives a correction of order N c when combined 
with the factor 1/f 2 . Thus it is not consistent with the 1/N C expansion. However this 
term has a trivial operator structure and therefore does not contribute to the baryon mass 
differences. Furthermore, because of mj, it is of 0(e 2 ) and suppressed in comparison with 
the leading order of the tree level mass. Secondly, the leading orders of O w , O k and O 71 are, 
respectively, 0(N 2 ), O(N^) and 0(N®). The leading order in 1/N C of each term is given in 
Table |T1[ One would expect that the gh and h 2 terms are suppressed as compared to the 



g term. This is true for the pion and kaon loop corrections as can be seen from Table pi 



However, in the case of //-meson loop, the gh and h 2 terms are as the same order as the g 2 
term. This is similar to what we have seen in the a 8 calculation in Section III. Thus in order 
to get the correct result for the rj loop corrections, we have to consider n-body operators in 
general, unless the coupling constants of such operators are numerically suppressed. In our 
estimate, we keep terms up to the 2-body operator in X m , i.e., the h term. Finally, we note 
that the g 2 terms involve the operators, {S}, {S}{S}, {cr J }{o r ' ? } and {«S<7 3 '}{cr- ? '}, which have 
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already appeared in the mass formula ( f4.18| ). This means that the g 2 terms satisfy the three 
mass relations, (Ml), (M2'), and (M5").0 Corrections to the mass relations come from the 
gh and h 2 terms which include {iSjlo"- 7 }!^}, etc. 

To estimate the loop correction from Fig. |4](a), we include the mass insertions to the 
intermediate baryon states. Let the mass difference be denoted by 8Mb'. Then the baryon 
self energy from this diagram reads 

E(u;) = -^(B\X^\B')(B'\X» a \B)l, u (uj) } (5.4) 

where 

J (27r) 4 \k ■ v — to + oM B ' J rn z aa — k z 
Calculation of the loop integral gives 

8M B = J 2 {5M B/ ,m u )^{B), (5.6a) 

where 

7 n (#) = Y,{B\X ia \B')(B'\X ia \B), (5.6b) 

a 

with a = 1, 2, 3 for the pion loop, a = 4, ... ,7 for the kaon loop and a = 8 for the eta loop, 
and 



X 3 



1 — In ( — J \, for m 2 A > x 2 , 



24vr 2 / 2 I V A J 



xm\ { , (m A \ 2 \ 1 , 9 9N ct /9l x-Jx 2 -m 2 A 
A J 2-31n — - \ + — (x 2 - m 2 A 3/2 In ■ V 



48tt 2 / 2 { V A / J 24vrV 2V " y x + y^T^ 

+ ^|l-ln(^) 2 }, for m\ <x 2 . (5.6c) 



In the limit 8Mb' = 0, we can recover the result ( |5.2|) . In the case of 8Mb' = (or constant), 
the loop correction can be represented in terms of the operators given in Eq. ( |5.2b| ). This 
is possible because the loop integral does not depend on the intermediate baryon state. 
However, this is not the case in Eq. ( |5.6|) since the loop integral depends on 8Mb 1 - 

We can write Eq. ( |5.6|) in a more convenient form as follows. With the usual definitions, 

a ±1 = T^(a x ±zay), a° = a z , (5.7) 



3 Note however that the mass relations (M2) and (M3) receive corrections from the g 2 term. 
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we use the Wigner-Eckart theorem, 

(V, f, m! |X(fc,g)| 7 , j, m) = ( _ 1)2 Jj^Mjj^O (y/||Z(fe)H T i). (5.8) 
Then after some algebra, one can show that 

Z a=l±i2,3 

7f(5) = f E [(s'll* a l|£)] 2 , 

Z a=4±i5,6±i7 

7 ^(B) = c B [(i?'||X 8 ||i?)] 2 , (5.9) 
where eg = 1/2 for octet baryons and 1/4 for decuplet baryons. Since 

x i+i2,i+i2 _ -gy/^alo^ - Ay2(QtJot 3 + a^a 4 )(ala 2 + ala 4 + c4a 6 ), (5.10) 

and so on, one can compute the matrix elements ^,(B) using the baryon wave functions 
given in Appendix A. The final results for ^,(B) are given in Appendix C. 
By comparison with Eq. fl5.2|) , we therefore have the relation 



O n (B) = l^B'(B), (5.11) 

which can be obtained by taking 8Mb 1 = in Eq. Q5.6Q . However, by inserting ^,(B) given 
in Appendix C, one can find that the above closure relation does not hold with B' 6 {8} 
and {10} only. This is because we have 

iy E \ B ')( B 'i (5-i2) 

B'={8},{10} 

in the large N c limit. The equality in the closure relation holds only for N c = 3. To form a 
complete set, we need an infinite number of states for infinite N c . However, fortunately in 
our case, X ta is a spin-1 operator. So what we need in order to satisfy the relation (|5.11|) 
is to include the intermediate baryon states up to spin 5/2. This is done in Appendix A, 
where we give all the states B' of spin 1/2, 3/2, and 5/2 to fulfill Eq. (|5.11|) . All these 
additional states are fictitious, i.e., they do not exist in the real world with N c = 3, but they 
are needed to satisfy the closure relation in the large N c limit. Note also that the baryon 
self-energy of Eq. (|5.6| ) starts at 0(e 2 ). 

The contribution to the baryon self energy from Fig. §(b) vanishes for the meson-baryon 
couplings ( |3.1D . The contribution of such a diagram comes from the effective Lagrangian 
( |4.2| ). Consider for example the one- loop correction from C^g of Eq. ( |4.2b| ) to the baryon 
self energy. This one-loop correction comes from the {m} term of C^, which is expanded 
as 

m = m i ~ 2^[ n ' [ n ' m J+]+ + • • • ' ( 5 - 13 ) 
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where [A, B], = AB + BA. Then the one- loop correction to the baryon self-energy reads 



2f 



r {[n, [n,m,] + ] + }A n , 



where 



A n = -,/ 



d 4 A; 1 



(2n) 4 m 2 u -k 2 ' 

By evaluating the loop integral using dimensional regularization, we get 



5M f 



where 



16vr 2 / 2 A 2 



(B\V?\B) = -^Y,{ll Xa ^™ q ) + ] + }. 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



In the same way, we can compute the baryon self energy of Fig. |](b) from the higher order 
terms of Eq. (f4.2|) to obtain 



6M. 



9 2 

mi, . m 



B 



S-ln^(5|P n |5), 



16tt 2 / 2 A 2 



(5.18) 



where 



|A a , [iA a ,mJ + ] + ) 



A-{[A a , [A<\m 9 ] + ] + a<}K} - | {mJA<\ [A°,m g ] + ] + } 
f ' 2 {mJ{[A a ,[A<\m g ] + ] + },-^tr 



47V t 

-^tr (m 9 [A<\[A a ,m 9 ] + ] 



^iV 2 



([A a , [\ a 



,™>q\ 



(5.19) 



Explicit calculation gives 



V" = --6 1 (2 J B m)[iV c - {S}} - 3(2B m) ai 

(2B m) Ua'Uo*} - {Sa l }{a 1 }} b 2 - 3(2B m) 2 [N c - {S}} c x 



2N r 



— (2B m)(2B ) [mN c + (m s - m){S}} [N c - {S}] c 2 
JL(2B )m{aJ}{a j }a 2 - 6(2B m) 2 p 1 , 



(5.20a) 



V K = --6i(2S )(m + m s )[N c + {S}] - 2(2B )(m + m s )a x 
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2N, 



{2B )(m + m s ) {a l }{a 1 } + {Sa l }{a 1 } 



-po 



- (2B ){m + m s ){2B Q ) [mN c + (2m s - m){S}] c x 

- —(2B )(m + m s ) [mN c + (m s - m){S}] [N c + {S}] c 2 

- ^{2B ){m s + m){a j }{a j }a 2 - 2(25 ) 2 (m s + m) 2 /3 1; 

--b 1 (2B )[mN c + (4m s - m){S}} - \{2B ){m + 2m s )a 1 
o 6 



(5.20b) 



1 



6N, 



(2B ) m{<x*}{(x 4 } + (4m s - m)^}^} 



- l -(2B,f [m 2 N c + (4m 2 s -m 2 ){S} 



Cl 



3N C 
1 

3N 2 



(2B ) 2 [mN c + (m s - m){S}] [mN c + (4m s - m){S}] c 2 
{2B )(2m s + rh){o 3 }{a j }a 2 - ^(25 ) 2 (2m 2 + m 2 )^. 



(5.20c) 



Thus the leading order of this loop correction is 0(e 4 ). 

However, there can be other one-loop corrections at 0(e A ) from higher order terms in 
the chiral Lagrangian, which can be written as 



5£ cS = A 1 {A»A IM } + ^{A,}{A»}. 



(5.21) 



Generally, terms which involve {(v ■ A) 2 } and {v ■ A}{v ■ A} are also possible. However, 
these terms can be absorbed into Eq. (|5.21| ) because of the following identity in dimensional 
regularization |i"4fl : 



d d k (vkf 



d d k k 2 



{2K) d m 2 -k 2 dJ (2vr) d m 2 - k 2 ' 
The Lagrangian ( |5.21|) gives the one-loop correction of the type of Fig. |](b) as 



SM 



B 



where 



Q n = E 



The leading order of this term is 0(e 4 ) since 



Q 
Q 



3A 1 



4 



(N c - {S}) + 



AN, 

2^ + ^ + M 



{a l }{a 1 } - 2{Sa i }{a i } + 2{S} + {S}{S} 
N c + (N c - 1){«S} + {Sa l }{a 1 } - 2{S}{S} 



Q v = ^(N c + 3{S}) + -^[N C - 3{S}] 2 . 



(5.22) 

(5.23) 
(5.24) 

(5.25a) 
(5.25b) 
(5.25c) 
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From the expressions for the operators V n and Q n of Eqs. (^20|) and QOBD , we can 
see that these are linear combinations of the operators that appeared already in Eq. Q4.18|) . 
This means that the loop corrections of Fig. [|(b) satisfy the mass relations, (Ml), (M2') 
and (MS'). 

In addition to the one-loop corrections of the previous calculation, we have to consider 
one more contribution, i.e., the 1 /Mb corrections. They have been calculated in Refs. 
H3yH|] within the framework of baryon chiral perturbation theory. To estimate the 1 /Mb 



corrections, one can use the relativistic form of the effective Lagrangian and then expand it 
to obtain the 1 /Mb terms. Or one may write down all possible next order terms in 1 /Mb 
and then fix the coefficients by using the so-called "velocity reparameterization invariance" 
The two methods should give the same result. 
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In this paper, therefore, we use the 
results of Ref. [2J| as discussed in Ref. [27] for a simple estimate on the 1/Mb corrections.^ 
If we consider the one-loop self energy of Fig. ^(a) with the intermediate state baryon mass 
Mb' in a fully relativistic theory according to Ref. P5fl , then we have 



6M 



B 



i(3 

w 



d*k j^/KI+I + Mb^ 
(2tt) 4 (k 2 - m 2 u )(2P ■ k + k 2 )' 



(5.26) 



where (3 stands for an SU(3) Clebsch-Gordan coefficient. By expanding the loop integral, 
one would have 



5M, 







B 



16np 



M%, 



7T 



7s + ln(4vr) + 1 - lnM|, 



+ n - - 1e + 1h(4tt) 

7r \ e 



2-lnM^, 



7lM BI 







--)] 


L . 





(5.27) 



where e = d — 4 in dimensional regularization. The first two terms proportional to M B , and 
Mb'Vp^ are the troublesome terms as noted by Ref. [ 25] . The m\ term is what we obtained 
previously, and the m^/MB' term is the 1/M B correction we want. Here we note that the 
1 /M B correction terms carry the same Clebsch-Gordan coefficient as the term. This 
was verified by explicit computation in Ref. [13 1. We use this result for our estimate of the 
l/M B corrections, 



5M 



B 




16n 2 f 2 M B 



(B\O u \B) 



(5.28) 



with O u defined in ( |5.2b| ). We can use M B = aoN c and note that the order of this 8Mb is 

0(e% 

Finally, we get the full one-loop correction to the baryon mass as 



4 See also Ref. [E26| for a critical review. 
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7D ^ f I 1 11 ^ 1 

5M B = J2j 2 (5M B ,,m u )(B\^\B) - |l + - ln^J (£|0 n |E) 



2 2 4 2 



where the operators, O u , V n and <2 n are respectively given in Eqs. (|5.3|), ( |5.20| ) and ( |5.25| ) 



and is defined in Eq. ( |5.6c|) . Note that when we calculate the 7^ term, we should 
include the fictitious intermediate baryon states of spin up to 5/2. From the structure of the 
operators, we can see that the loop corrections to the mass relations (Ml), (M2') and {MS') 
come from the 7^/ and 1/Mb terms, and the other terms respect the three mass relations. 
Note also that the leading contribution to 7§/ is 0(e 2 ) while those of O n , V n and Q n are 
0(e*). 



VI. SIGMA TERM AND STRANGENESS CONTRIBUTION TO THE NUCLEON 

MASS 

The pion-nucleon sigma term, defined as 

oVtv = m(p\uu + dd\p), (6.1) 

can be computed from the expression of the nucleon mass using the Feynman-Hellman 
theorem: 

- 8Mn (r o\ 

a nN = m— — . (6.2) 

am 

The strange quark contribution to the nucleon mass can be written as 

(p\m s ss\p) = m s -7^- ■ (6.3) 

Then we can estimate the strange quark matrix element (SME) (p\m s ss\p) from the mass 
formulas derived in the previous Sections. 

In this Section, we consider the SME at the tree level. Up to O^e 1 ), the nucleon mass is 
written as 

M N = a N c + —ai + Nju 2 ^. (6.4) 



We find that there is no strange quark contribution to the nucleon mass at this order: 

v-kN = Ncmlh, 

(p\m s ss\p) = 0. (6.5) 

From Table [111 we observe 

a = a N c + N a mlbx = 1063 MeV, b = ^ = 26.2 MeV, (6.6) 



N c 
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where a and b are defined in Eq. (|4.18|) . So using a n N = 45 MeV f28f , we can fix the three 
parameters as 

a = 339.3 [337.7] MeV, a x = 78.6 MeV, 6j = 7.88 [8.75] x 10" 4 MeV -1 , (6.7) 

where the values in square brackets correspond to a n N = 50 MeV as suggested by the lattice 
calculation of Ref. ||2"9"|j . 

The non- vanishing SME comes from the 0(e 2 ) terms. The nucleon mass up to this order 
reads 

M N = a N c + — a x + Njn 2 ^ + (2m 2 K + m\)a x , (6.8) 

and involves four parameters. We find 

o-kN = ml(N c bi + 2qi), 
(p\m s ss\p) = (2m 2 K - m^cm, (6.9) 



which gives 

(p\m s ss\p) = l -(2m 2 K - ml) I ^ - \, l>i ) . ( O.J 0) 



Note that the SME starts at 0(N®) in 1/N C counting as pointed out in Ref. |T0[ . From the 
best fit of Table 0, we get 



a = a N c + Ncmlh + (2m 2 K + m\)a x = 923.9 MeV, 
b = — = 19.54 MeV, 

c = 2(m 2 K — m 2 )bi = 159 MeV, (6.11) 



which gives 



a = 190.45 [168.05] MeV, a x = 58.62 MeV, 

b x = 3.52 x 10~ 4 MeV" 1 , a x = 6.53 [7.85] x 10" 4 MeV" 1 , (6.12) 

for a n N — 45 MeV (the values in the square brackets are for a^N = 50 MeV). Then we have 

(p\m s ss\p) = 307.8 [369.6] MeV. (6.13) 

This shows the familiar strong dependence of (p\m s ss\p) on the value of u^n- This is 
because the constant multiplying ova? in Eq. ( |6.10|) is as large as 12.4. For example, if we 
use a nN = 65 MeV, we find 555 MeV for the SME. 

However, we have to include at least the 0(e 3 ) terms to get a more reliable value of SME 
because the fitted baryon mass spectra is reasonably consistent with the experiment from 
this order onward. For the nucleon mass we have two additional terms so that 

3 3 
M N = a N c + — oi + N c mlh + (2m 2 K + m^ai + j^m% + m 4 A^ c ( Cl + c 2 ). (6.14) 
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Although there are altogether 7 parameters in the Lagrangian, we have only 6 independent 
parameters since c\ and c 2 enter in the form {c\ + c 2 ) for all baryon masses. The final result 
is: 

3 

a vN = ml[N c b x + 2ai + — b 2 + 2iV c m*(ci + c 2 )], 
(p\m s ss\p) = (2m 2 K - m\)a x , (6.15) 

which implies 

(p\m s ss\p) = l -(2m 2 K - m£) j^J - - J-6 2 - 2iV c m£(ci + c 2 )| . (6.16) 

To estimate this matrix element, we must determine the parameters. Not all of them can 
be fixed, however, since there are 6 parameters while we have only 5 pieces of information: 
four from baryon masses and one from the irN sigma term. From Table |T1], we have 

a = a N c + Njnlbi + {2m 2 K + m*)ai + A c m 4 (ci + c 2 ) = 863.7 MeV, 
2 

b = W + jt h = 25 '° MeV ' 



c = 2(tt4 - m l) b i = 227 - 8 MeV , 
2 



d = ^-{m 2 K - ml)b 2 = -16.6 MeV, (6.17) 



which gives 

ai = 77.03 MeV, 61 = 5.04 x 10" 4 MeV" 1 , b 2 = —1.10 x 10 -4 MeV -1 . (6.18) 

Note that these best fit values of a± and b± at 0(e 3 ) are between the values found at 0(s 1 ) 
and at 0(e 2 ). Since the other parameters cannot be determined uniquely, we rewrite the 
SME of Eq. (|1D in the form: 



where we have expressed (ci + c 2 ) in terms of a^N and a. Since a is fixed by the mass 
spectrum, therefore, the SME of the above form depends on the unfixed parameter clq. For 
a numerical estimate we can use the fitted values of ao from the calculations at 0(e 1 ) and 
0(e 2 ), i.e., ao = 190 ~ 340 MeV. This leads to {p\m s ss\p) ranging between about 250 MeV 
and —190 MeV. Now the dependence on the ttN sigma term is very weak, while it depends 
strongly on the value of ao, leaving {p\m s ss\p) almost completely uncertain. 

At 0(e 4 ) and 0(e 5 ), the situation becomes even more subtle. There are 9 parameters 
with 5 pieces of information in case of 0(e 4 ). If we take into account the corrections from 
0(e 5 ), then we have 13 effective parameters^ with 7 constraints. Additional information is 



5 There are totally 15 parameters up to 0(e 5 ) calculation. However, the three parameters di, 2 ,3 
appear only in the form of (d\ + d 2 + ^3) at this order. So there are effectively 13 parameters. 
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therefore required such as isospin symmetry breaking effects in the baryon spectra and/or 
KN sigma terms ||30|| . As a reference, we give the formulas of the sigma term and the SME 
up to 0(e 4 ) below: 

<t*n = ml |iV c 6i + 2«x + J-6 2 + 2N c ml( Cl + c 2 ) + -^a 2 + 4m*/3i j , 

(p\m s ss\p) = (2m 2 K - m 2 n ) L, + j^a 2 + 2{2m 2 K - rr%)h\ , (6.20) 



where 



a N c + Njnlh + (2m 2 K + m 2 )^ + N c m*{c x + c 2 ) 
+ (Am 2 K - Am 2 K ml + 3m*) fa = 863.7 MeV, 

aT + Tfh + TH( 2m K + ml)a 2 = 25.0 MeV, 



iV c N c iV c 2 
c = 2(m 2 K - m%)b x + 4(m^ - m\)[m 2 K c x + m^ca] = 227.8 MeV, 

d= ^r{m 2 K -ml)b 2 = -16.6 MeV. (6.21) 

In essence one observes that corrections beyond the standard estimate ( 6.1 0|) for (p\m s ss\p) 
are so large that they prohibit quantitative conclusions about the strange quark contribution 
to the nucleon mass. 



VII. SUMMARY AND DISCUSSION 

In summary, we have re-analyzed baryon masses within baryon chiral perturbation theory 
in combination with the large N c expansion. Before computing the baryon masses, we have 
calculated the baryon axial current. We find that the two diagrams of Fig. [TJ give contri- 
butions of the same order in 1/N C counting. Inclusion of the wave function renormalization 
terms is crucial to get the right order for the one-loop corrections because this gives the 
proper commutator structure to the baryon axial current operator. However, when calculat- 
ing a%t B , two-body operators give contributions of the same order as one-body operators. 
Unless the coupling constants of the general n-body operators are suppressed numerically, 
their effects must be included in order to be consistent with the 1/N C expansion. 

Next, we have considered the baryon mass spectrum in this scheme. For this aim, we 
have used that both mn and 1/5M scale as 0(e), where mn and 5M, respectively, represent 
the Goldstone boson mass and the octet-decuplet mass splitting which depends on l/N c . At 
the tree level, we found that the empirical mass spectrum is well reproduced to 0(e 3 ) and 
the corrections from 0(e 5 ) are not so crucial. But the Gell-Mann - Okubo mass relation and 
the equal spacing rule in the decuplet are modified at 0(e 5 ). At the one-loop level, there is 
no additional contribution that gives the characteristic commutator structure, and the loop 
corrections seem to violate the 1/N C expansion. However, the leading terms are constant 
for all baryon states and can be safely absorbed into the central baryon mass in the chiral 
limit. The meson loop corrections involving the operators O u with the coupling constant g, 
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V u and Q n of Eq. (]Q9|) satisfy the modified mass relations (Ml), (M2') and (AfS')- To 
get the correct result, the intermediate baryon states must include fictitious states of spin 
up to 5/2 in order to satisfy the closure relation, J2b \B)(B\ = 1; f° r the spin-1 operator 
X m in the large N c limit. As in the calculation of a%> B , the r^-meson loop corrections to the 
baryon self energy require general n-body operators in order to be consistent with the 1/N C 
expansion. 

Finally we have estimated the strangeness contribution to the nucleon mass at the tree 
level. We confirmed that this matrix element is 0(N®) in the 1/N C counting. At leading 
order, namely 0(e 2 ), this contribution can amount to more than 300 MeV. At the next 
order, we cannot uniquely determine the mass parameters because of lack of independent 
empirical information. But the upper bound of the strangeness contribution to the nucleon 
mass is now reduced to around 250 MeV. 
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APPENDIX A: BARYON STATES 



The octet and decuplet baryons states \B,j z ) in the number space are given in this 
Appendix. For simplicity we give only the s z = +1/2 states for baryon octet and s z = +3/2 
states for baryon decuplet. Other spin states can be obtained straightforwardly. The octet 
states are 



\p, 


+\ 


\n. 


+\ 


|A. 








|E°, 




£~ 






+1 







C N a\ (4)10), 
C N at(Al) n \0), 
-C A al (4)^0), 
-C S «I (4) (Atr^\0), 

-^C E {aI4 + «!4K4) 

c E 44 (4) n-1 |o), 
-c H 44(4) n_1 |o), 



(Al) 
(A2) 



n-1 



10), 



(A3) 
(A4) 



where 4 
constants 



a\a\ — a\al, A\ = a\a\ 



\n\C 



a\a\ and A\ = a\a\ — a^al, with the normalization 

1 



N 



[(n-l)!C E 



n + l)(n + 2) 
,2 2 



\n\C 



A I 



n + 1' 



n(n + l)(n + 2)' 



[(n-l)\C B ] 



(A5) 



3n(n + 1) ' 



23 



from the condition (B,j z \B,j z ) = 1, where N c = 2n + 1. The negative signs of some 



states were introduced to be consistent with the quark model convention j3T|. Explicitly the 
spin-up proton state can be written as 



|p,+§) = CV£ 



nl 



(-1) 



n 



k + 1, k, k, n — k, 0, 0), 



(A6) 



by making use of 



(a+bt=y: 



nl 



k=0 



k\{n-k)\ 



j^n—k j^k 



(A7) 



The decuplet states are as follows: 



|A++ 


+§) 


= C A a\a\ a\ (AD^O), 


|A + 


+1) 


= V3C A a\ a\ a\ {A\) n ~ x \Q\ 


|A°, 


+1) 


= V3C A 4 44 (AD^lO), 


|A~, 


+1) 


= C A ata&ai (4) n_1 |0), 




+1) 


= Cx* 444 (4)" _1 |0), 


|S*°, 


+1) 


= y/2Cv* 4 alaUAlY^lO) 


|E*~, 


+1) 


= 0^444 (4) n_1 |o), 




+1) 


= cu 444 (Air- l \o), 




+1) 


= cu 444 (Air- l \o), 




+1) 


= c* n 444 (4) n_1 |o), 



(A8) 



(A9) 

(A10) 
(All) 



where 

[{n-l)\C A } 2 - 

\{n-l)\C^] 2 



n(n + l)(n + 2)(n + 3)' 
1 



[(n-l)!C a .] 



n(n + l)(n + 2) 



n(n + 1) 



\(n-l)\Cn\ 2 = — . 



(A12) 



However, the octet and decuplet states are not sufficient to satisfy the closure relation 
( |5.11| ) in the large N c limit, and we have to include higher spin states. Since X ia is a spin- 
1 operator, it is sufficient to introduce fictitious states up to spin 5/2. These states are 
distinguished from the octet/decuplet by a tilde and we denote the strangeness S = —4 
states by |S). These states can be obtained by considering 5-quark states multiplied by 
(Al) n ~ 2 , whereas the conventional octet and decuplet are formed by 3-quark states with 
(Al) n ~ 1 . Then the fictitious states of spin 1/2 are 



F2, +o) 



ck«i44(4r 2 |o), 

i=Q {24 A{ + 4 A\}A\{A\r 



! |o), 
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-3- 



+h) = 



{44 + 244} 4 (4r 2 |0), 



| "4 

l«2,+i) 



|) = Q444(4r 2 |o) 



~) = 0,444(4r 2 |o), 
= v^c>i44(4r 2 |o), 
= Q.444(4r 2 |o), 



where 



[(n 2)!Cy 2 - („_!)„(„+!)(„ + 2) ' 

""- 2 ) !C "' 2 = 4(n-l)!,(n + l) - 
For the spin 3/2 states, we have 

|S;,+|)=C7 £ . 4444(4) n ~ 2 |0), 

\% +§) = \ C ±* {4 4 4 A{ + 3a\ 4 4 4} (4) n ~ 2 |0), 

ft +|) = ^fc7 £ . {4 aj 4 4 + «I 4 4 4} (4) n " 2 |0), 
-|) = \ C±* {&*} at at A{ + aj 4 4 4} (4) n " 2 |0), 

-D = c^4444(4r 2 |o), 

-D = c7g. 4444(4r 2 io), 

, +|) = -L <? g , {4 a} at 4 + 24 4 a\ 4} (4) n " 2 |0), 
+ |) = ^= <?§. {2al 4 4 4 + «! 4 4 4} (4) n ~ 2 |0), 

D = c7g. 4444(4r 2 |o), 
" =^4444(4r 2 |o), 



where 



i — * 

"2 



"3 

^* 
"4 

fit 



+ 
+ 

+ §) = ^= Qj. {«1 «5 4 4 + 4 4 



+|) = c7 fi . 44 4 4(4r 2 |o) 



^4 4} (4)' 



—2 1 



2/ a z o o o w \ s/ i / ■ 

|s;,+|) = c7 8 . 4444(4r 2 |o), 
fe+|) = c7 8 . 4444(4r 2 |o), 



2 4 

[(n-2)!C f: .] = ^— 1)n ( n + 1) ( n + 2 )( n + 3) 
[(n-2)!C § ,] = 5 (n-l)n(n+l)(n + 2)' 



(A13) 



(A14) 



(A15) 



(A16) 



(A17) 



(A18) 
(A19) 
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[{n-2)\C^f = 
[(n-2)!Cy 2 = 
The possible spin 5/2 states are 

iAr,+§) 



5(n — l)n(n + 1) ' 



I A ** 

I A ** 
1^3 

I A ** 
1^4 



I AS 



^3 
^4 
^5 

| 

I Z?** 
I "2 

I ^** 
I "3 

1 1— ■•** 

ins* 

isr 
isr 



+1) 
+1) 
+1) 
+1) 
1) 

2^ 
2^ 

2> 

^) 
2> 

-) 

2^ 



15(n — l)n 

,t „t „t „t „ ; t (4) 



C A « a{ a{ a\ a\ a\ (AJ) n -^|0), 



v 7 ^ C A » aj a} aj a} a\ (Aj) n ~ 2 |0), 

a/5 C A .. a} a} «}«!«! (4)»- a |0), 
VlOC A „aIaIaJaJaS (A\) n - 2 \0), 
V5C A ..aJaSaSaSaS(Al)"- 2 |0), 

c A « 44444(4r- 2 |o), 

a\ a\ a\ a\ a\ (4) n ~ 2 |0), 

2 (%„ aj a{ a\ a\ a\ (Al) n ' 2 \0), 
t ™t „t „,t „t 



a{ a| «3 4 (4) n " 2 |0), 



2 Cg» 



a| a: 3 a\ a\ a\ 



+§) = 
+ §) = 
+§) = 
+§) = 

+§) = 
+§) = 
+§) = 
+§) = 



(Aj)»- 2 |0), 

C S « 4 «3 «3 «3«5 (4)"" 2 |0), 

atajajaj (4) n ~ 2 |0), 

ra-2 



>/3 Cg.. a} a} a| ^ ^ ( A\)' 



|0), 



V3C§..a}aSa!«!«J(^) B " 2 |0), 

a«44444 (4) n_2 |o), 

C„„ aM 44«5 (4)"~ 2 |0), 
v^C fi ..a}aSaJaJaJ (At)- 2 |0), 



i(4) n - 2 io), 



CK3 «3 CK5 0:5 ct 



„t „t „t „t „t (4r 2 | ), 



a} 0:5 CKj a§ a§ 
at at at at at 



(4) n ' 2 |o), 



(A20) 



(A21) 



(A22) 



(A23) 

(A24) 
(A25) 



where 



[(n 


— 2)!C A **] 


[(n 


-2)!C fi ..] 2 


[(n 


- 2)!Cs„] 2 


[(n 


-2)!C7 fl ..] 2 


[(n 


— 2) !Cg»* ] 



(n - l)n(n + l)(n + 2)(n + 3)(n + 4) 
5 

(n- l)n(n + l)(n + 2)(n + 3) ' 
2 



(n- l)n(n + l)(n + 2)' 
1 

2(n- l)n(n + 1)' 
1 



12(n- l)n" 



(A26) 
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Note that the A, S, H, Q and S families have isospin 5/2, 2, 3/2, 1 and 1/2, respectively, 
and their normalization constants contain the factor l/(n — 1) so that these states can not 
be defined with N c = 3. Using the results given in Appendix C, one can find that these 
fictitious states ensure the relation ( |5.11|) . 



APPENDIX B: EXPLICIT RESULTS OF 



BB' 



In this Appendix, we give the explicit results of 0^, and P^, of (|3.22| ) from the g term 

of ]•:<)• » 



and 



ol+i2,n 
Ppn 

nl+i2,K 
Ppn 

ol+i2,r) 

r^pn 

Pay,- 



nl+i2,K 

Pat,- 



nl+i2,r] 

Pas- 



nl+i2,n 



nl+i2,K 

Pe°e- 



PS 



l+i2,ri 



0: 



l+j2,7T 

E°E" 



M+i2,K 

Pe°e- 



P 



l+«2,?7 
E°E" 



-l(N c + 2)g 3 -~(N c + 2)g, 
-l(N c + 2)g 3 -±{N c + 2)g, 



9 



(N c + 2)g 3 



3^2 
1 

2\?2 
1 

9\7I 
2N C 



(N c -l)(N c + 3)g 3 
(N c -l)(N c + 3)g 3 
(N c -l)(N c + 3)g 3 



1 



3V2 
1 

6Vl 



(A c -l)(A c + %, 
(A c -l)(A c + 3)<?, 



9~ 9 9 9 ' 



N 



.Ac s 
27 9 ' 



N r . 



;9, 



3^2 
1 

2~V2 
1 



(A c +iy 
(N c +l)g 3 
(N c + l)g 3 



1 



3^ 
1 

6\/2 



(A c + 1)^, 
(A c + 1)<7, 



(Bl) 



(B2) 



(B3) 



(B4) 



P 



4+i5,7r 



pA 



PpA 

Pae- 



Pit 



11 



isV A c + 3^ + - V / A c + 3^, 



9 



1 



16^3 v'^- 1/ -^K-l) 9 , 

o / o4+i5,7r 



(B5) 
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fig" = -^fV^V - yW-l)9, (B6) 



04+45,^ _ r>M+i5,7r 
-pSO ~~ Z ' J pA ' 

^ = -t^V^V - iv^T^ ( B7 ) 



144 v ° 16 
5^ /- , 5\/3 



M+ib,K _ M+i5,ir 



y pA 



55^3 r- , 5^ 



= -^V^ + ~^Nc + 35- (B8) 



For p 18 , we have 



7 l(iV) = 0. 



(B9) 



/3 8 >* = — o 3 8 8 ' K = — q 3 - — - — q 8 8 ^ = — q 3 

P PP 2 ^9 , P PP ^9 ^9, P PP Q ^9 , 

A? 8 ' 71 " - ^ 3 a8,A" _ 7 3 ^ 3 8 j) _ 2 3 

^ S " VS 9 ' ^ ~ 6VS 9 2y/3 9 ' ' ^ ~ 3VS 9 ' 

I&5 = ^9\ m = -^r9 3 + -^=9, = -^r9 3 . 

The constants are the same as the g 3 terms of 

APPENDIX C: MATRIX ELEMENTS OF j%,(B) 
In this Appendix we give the matrix elements of j^,(B). 

rtW = 4^^^ + 2 ^ + 3h ?> = \{N C - 1)(N C + 5)«? 2 , 

7f W = ^^[iV c , + /f , 7 f(iV) = ^[iV c ,-3/f , (ci) 

7 ^(iV) = (iV c -iy, T^^b + ^f, 
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7l(A) = 0, 

lUV = \(N c -l)(N c + 3)g 2 , 
7f(A) = 0, 
7 |(A) = 2(iV c -l)^, 
7§(A) = o 



1 *(A)= 1 -(N c -l)(N c + 3)g 2 , 
7£(A) = ^^[A^ + /f, 

7 g,(A) = 0. C 



7l(X) = -(iV c -l)(iV c + 3)£ 2 , 

7£.(S) = ^(iV c +l)V, 
7^(S) = ^(iV c + % 2 , 
7i(S) = ^ c + % 2 , 



7^(S) 



4JV? 



[27V c <7 + (7V c -3)/i] 2 , 



7S(S) 



7£(£) = ^[A^-3/f, 

7 " (s) = ^f [5iv ^ +3/l]2 ' 

tX(s) = o, 

7M£)=2^ 2 . 



[7V 2 5 - 9/f , 



[A^ + 3/f, 



7f (S) = ^ 
7 f«(S) = i(iV c + 3)/, 
im = ^N c g-{N c -Q)h}\ 



7S.(S) = ^ 2 , 

7^) = ^[5iV c , + 3/f, 

7^(S) = (iV c + lk 2 , 
7l*(S) = 2^ 2 . 



7 -(A) = -(AT c -l)(7V c + 5)<? 2 , 
7f (A) = 0, 

7 £(A) = ^^[A^ + 3/f, 
7l(A) = ^ + /f- 



7l(A) = i ^[iV c (iV c + 2)^ + 15/ i ] 2 , 
7 f(A) = ^ c + % 2 , 
7^(A) = 0, 
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7l(S*) = 1(JV C - l)(iV c + 3)/, 7 S(S*) = ^(iV c + 1) V, 

7£.(S*) = ^ W^c + + 12/.] 2 , 7^(S*) = \(N C - l)g\ 
llV*) = 5 -%^i N c9 + 3/f , tT(S-) = i(* c + % 2 , 

7i(S*) = ^^[^<7 + 3/f , 7^*) = 0, 

z ±iV c 

7S(S*) = \Nt g \ 7g.(S*) = ^r[^ + 9/f , 

7f (3*) = ^(^Vc- 1)<? 2 , t£(3*) = i(iV c + 3)^ 2 , 

7&(3') = W^[N e g + Sh]*, 7^*) = ^j^[N c9 + m\ 

72(S*) = <? 2 , 7 2*(H*) = -£=[N e g - (N c - 6)/f . 

7n(«) = 0, 7^) = (iV c + 1)<? 2 , 

= 5( 4^V cg + 3/f, 7 S(«) = ^[SAT^- (iV c - 9)/f . 
For the fictitious intermediate states, we have 

7i.(S) = t|(^ " 3)(iV c + % 2 , 7 f (E) = 2( ^ 3 W - 3/f, 
7 f,(S) = y(iV c -3), 2 . 

7l(S) = - 3)(iV c + % 2 , 7 |.(S) = A (JVc - 3)(iV c + % 2 , 



£(3) = { -%^lN c9 + 3hf, 7 ^(S) = f (iV c -% 2 . 



c 



7i..(A) = ^(iV c - 3)(iV c + l)g\ 7 ^(A) = 3( ^ 3) [^g - 5/f , 
7^(A) = ^(iV c -% 2 . 

7 |.(E*) = ^(N c - 3)(N C + % 2 , 7 5..(E*) = jj(JV c - 3)(iV c + 5)/, 

Tf(S') = y^(iV c - 3)^ 2 , 7 £(E*) = ^[iV cS - 15/f , 
7^(S*) = ^(iV c -3)^ 2 . 



c 
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7i(S*) = l(iV c - 3)(7V C + 3)/, 7l*(S*) = ^ " 3)(W C + 3)/, 

7l»(S*) = jjW* - 3)(iV c + 3)g\ 7 ^(S*) = 1(N C - 3)g\ (C13) 

7£(S*) = ^f-l N *9 + 15^] 2 , 7&.(3') = l(N c - 3)/. 



7 §(n) = ^(N c - 3)(N C + 1)/, 7 g.(n) = ^(iV c - 3)(7V C + 1)/, 

7§«(n) = ^(iV c -3)(iV c + l), 2 , ^(fl) = + , (CM) 



c 



7 s ^) = 3(7V c -3)(? , 

and the others are zero. Note that all the matrix elements with fictitious intermediate state 
contain the factor (N c — 3) so that they vanish in the real world with N c = 3. 
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FIG. 1. One-loop corrections to the baryon axial current. 
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FIG. 2. Wave function renormalization of one-loop. 
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FIG. 3. Best fit of baryon masses (tree) up to 0(e 5 ). Thick lines represent degenerate states. 




(a) 



(b) 



FIG. 4. One- loop corrections to the baryon mass. The filled-triangle denotes the mass insertion 
to the intermediate baryon state and the filled-box represents the meson-meson-baryon-baryon 
coupling from the chiral Lagrangian of ( f4,2|) and ( 5.21| ). 
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TABLES 





{U} 




{S} 




{Saj}{aj} 


{Ua 3 }{Sa 3 } 


N 


2n + l 


3 














A 


2n 


3 


1 


1 


3 





£ 


2n 


3 


1 


1 


-1 


-4 




2n- 1 


3 


2 


4 


4 


-4 


A 


2n+ 1 


15 














S* 


2n 


15 


1 


1 


5 


2 


"* 


2n- 1 


15 


2 


4 


10 


2 




2n - 2 


15 


3 


9 


15 





s 


2n- 1 


3 


2 


4 


-2 


-10 




2n- 2 


3 


3 


9 


5 


-10 


s* 


2n 


15 


1 


1 


-3 


-6 




2n- 1 


15 


2 


4 


4 


-4 


ft* 


2n - 2 


15 


3 


9 


11 


-4 


s* 


2n- 3 


15 


4 


16 


18 


-6 




2n+ 1 


35 
















2n 


35 


1 


1 


7 


4 




2n- 1 


35 


2 


4 


14 


6 


ft** 


2n - 2 


35 


3 


9 


21 


6 


g** 


2n- 3 


35 


4 


16 


28 


4 



TABLE I. Matrix elements of various operators for baryon states 
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Particle 


0{e) 


0(e) 


U(£ ' ) 


0(e°) 


Expt. 


N 


1142 


982 


939 


937 


939 


A 


1142 


1141 


1117 


1119 


1116 


Zj 


1 1 A O 

1142 


1141 


1183 


1183 


1193 




il4z 


1 Qnn 
ioUU 


lozo 


loz / 




A 

A 


1456 


1217 


1238 


1236 


1232 




1456 


1376 


1383 


1386 


1385 


i— <# 


1456 


1 CT O CT 

1535 


1528 


1530 


i cr or» 

153U 




1456 


1694 


1673 


1570 


1672 




424 


79 


16 


15 




a 


1063.0 


923.9 


863.7 


862.4 




6 


26.2 


19.5 


25.0 


24.9 




c 




159.0 


227.8 


96.5 




d 






-16.6 


51.8 




e 








-70.4 




f 








-67.8 





TABLE II. Best fit of baryon masses (tree) in the unit of MeV at each order of e using the 
formula (|4.18j ). 



operator 


g 2 term 


g/j term 


h 2 term 




N 2 


A° 


iV 6 7 2 


Q K 


A 1 








A 

c 


NO 


ATO 



TABLE III. The leading order of operator O depending on the coupling constants. 
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